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ABSTRACT 



X-rays from planetary nebulae (PNs) are believed to originate from a shock driven into the fast stellar- 
wind (v ~ 1000 km s ) as it collides with an earlier circumstellar slow wind (v ~ 10 km s ). In theory, 
the shocked fast wind (hot bubble) and the ambient cold nebula can remain separated by magnetic fields 
along a surface referred to as the contact discontinuity (CD) that inhibits diffusion and heat conduction. 
The CD region is extremely difficult to probe directly owing to its small size and faint emission. This has 
largely left the study of CDs, stellar-shocks, and the associated microphysics in the realm of theory. This 
paper presents spectroscopic evidence of ions from the hot bubble (kT » 100 eV) crossing the CD and 
penetrating the cold nebular gas (kT m 1 eV). Specifically, a narrow radiative recombination continuum 
(RRC) emission feature is identified in the high resolution X-ray spectrum of the PN BD+30°3639 indi- 
cating bare C VII ions recombine with cool electrons at kT e = 1 .7 ± 1 .3 eV. An upper limit to the flux of 
the narrow RRC of H-like C VI is obtained as well. The RRCs are interpreted as due to C ions from the 
hot bubble of BD+30°3639 crossing the CD into the cold nebula, where they ultimately recombine with 
its cool electrons. The RRC flux ratio of C VII to C VI constrains the temperature jump across the CD 
to AkT > 80 eV, providing for the first time direct evidence of the stark temperature disparity between 
the two sides of an astrophysical CD, and constraining the role of magnetic fields and heat conduction 
accordingly. Two colliding-wind binaries are noted to have similar RRCs suggesting a temperature jump 
and CD crossing by ions may be a common feature of stellar-wind shocks. 



Subject headings: planetary nebulae:individual: BD+30°3639 
flows - stars: Wolf-Rayet - X-rays: stars 
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1. Introduction 

High velocity interactions among different as- 
trophysical media due to collective plasma effects 
on length scales much smaller than the mean free 
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paths for particle collisions are commonly referred 
to a s collisionless s h ocks. For extensive reviews 
see iTidman & Kralll dl97ll) : iMcKee & Hollenbach 
These shocks have been suggested to account 
for a wide variety of astrophysical phenomena ranging 
from the Earth's weak bow shock to strong shocks in 
supernova remnants (SNRs). Despite elaborate the- 
oretical investigations, the detailed physics of colli- 
sionless shocks is far from being understood, mostly 
because direct observations or physical measurements 
of these shocks are notoriously difficult due to the ex- 
tremely small size scales on which they occur. 
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The extended X-ray emission in planetary nebu- 
lae (PNs) is most commonly explained by a (reverse) 
shock driven into the fast wind expelled by the cen- 
tral star as it collides with the earlier slow and mas- 
sive wind (see review by Frank Il999h . High angu- 
lar resolution X-ray images of PNs seem to support 
this notion as the X-ray gas referred to generally as the 
hot bubble appears confi ned to the interior o f the opti- 
cal and IR nebula (e.g.. lKastner et al . 2000). The low 
X-ray temperatures and observed luminosities cannot 
be explained by the present day wind velocities and 
mass loss rates. This perhaps suggests that the X-rays 
are due to the fast wind ejected during the early-PN 
phase or late post-asymptotic giant branch (post-AGB) 
phase, when the rapidly evolving wind was slower 
(v ~ 500 km s" 1 ) and its mass loss rate was higher than 
what it is today dAkashi et alj 12006). Alternatively, 
conduction of heat from the fast wind to the slow wind 
has been invo ked to explain the observed lo w X-ray 
temperatures dSokerill994llSteffen et ai1l2008l) . In that 
case, it would be the nebular gas evaporating into the 
hot bubble that emits the X-rays. The fast PN wind 
running into the denser and slower nebular gas leads 
to relatively slow shocks (v, * 30 km s _1 ) propagat- 
ing in this dense gas. In a few cases, the morphology 
of the PNs sug gest that the fast wind is at least m ildly 
collimated dKastner et alj|2003t ISahai et alj2003b . 

High resolution X-ray spectroscopy of PNs is par- 
ticularly challenging owing to the low X-ray flux from 
these sources. The only published X-ray grating ob- 
servation of a PN is the ChandrafLETG S 300 ks ex- 
posure ofBD+3023639 dYu et alJl2009albT) . The anal- 
ysis of lYu et al.1 d2009al) reveals a temperature range of 
(at least) 150-250 eV and extremely non-solar abun- 
dances. The elements C and Ne appear to be particu- 
larly enriched with respect to solar abundances (C/O ~ 
30, Ne/O ~ 4), while Fe and N are deficient (Fe/O ~ 
0.2, N/O ~ 0.4). In fact, these abundances found in the 
X-ray plasma closely trace the abundances measured 
directly from the wind of the present-day WC central 
star of BD+30°3639 indicating that the X-rays do in- 
deed originate from the fast wind, while evaporation of 
nebular-composition gas and possibly also conduction 
of hea t to the surrou nding nebula have no observable 
effect dYu et al.ll2009al) . 

The present paper is dedicated to an unusual emis- 
sion feature in the LETGS spectrum of BD+30°3639 
around 25.30 A, which we interpret as a narrow ra- 
diative recombination continuum (RRC) of bare C VII 
forming H-like C VI. This feature, observable only at 
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Fig. 1 . — LETGS background-subtracted spectrum of 
BD+30°3639, using bin size of 0.1 A, overlaid with 
the best-fitted model 



the high spectral resolution of the LETGS, provides a 
new insight into the region of the shocked fast wind of 
BD+30°3639 and the microphysics at play. 

2. Observations and Model 
2.1. The Spectrum 

BD+30°3639 was observed with Chandra for a to- 
tal of ~ 300 ks, broken into several segments from 
February - December 2006. All observations were car- 
ried out with the low energy transmission grating spec- 
trometer (LETGS) and advanced CCD imaging spec- 
trometer (ACIS-S) configuration. The data were pro- 
cessed using the standard Chandra pipeline and com- 
bined into one spectrum using the standard tools of 
the CIAO software package. The full details of the 



obser vations and data reduction are given in lYu et al. 
d2009ah . 



As an initial step and following lYu et al.l (|2009a), 
we employed the XSPEC package (version 12.3.1, 
lArnaudlll996l) to fit the entire spectrum with two tem- 
perature components (2 -T) of the astrop hysical plasma 
emission code APEC dSmith et al. 2001 ). photoelec- 
trically a bsorbed by cold gas usi n g the absorption 
model of Morrison & M cCammon J 19831) . A minor 
improvement over the Yu et al.l \ 2009ah approach is 
introduced for modeling the spectral line profiles; 
wavelength-dependent Gaussian smoothing is used. 
Since BD+30°3639 is ~ 4" across and spatially re- 
solved by Chandra, and since LETGS is a dispersive 
slitless spectrometer, emission lines appear broadened 
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Table 1 : Parameters of the best-fitted model 



Parameter 


Value ± Error 


Yu et al. (2009a) 1 ' 


Units 


N* 


0.24 ± 0.04 


0.24 ± 0.04 


10 25 cm" 2 


Thermal Components 


kT\ 


260 ± 30 


250 ± 30 


eV 


EM l 


50 ± 10 


70 ±50 


10 53 cm" 3 


kT 2 


160 ± 30 


150 ± 30 


eV 


EM 2 


160 ± 45 


170 ± 35 


10 53 cm" 3 


Ac/Ao 


37 ±8 


33 ± 15 


solar 1 ' 


An/Ao 


0.3 ± 0.4 


0.4 ± 0.4 


solar 1 


Ac, 


0.9 ± 0.4 


0.9 ± 0.4 


solar 1 


A N JA 


3.8 ± 0.4 


3.8 ±0.8 


solar 1 


AuglAo 


0.8 ± 0.5 


0.7 ± 0.5 


solar 1 


Asi/Ao 


2.0 ± 1.1 




solar 1 ' 


A F JA 


0.15 + 0.1 


0.2 ± 0.2 


solar 1 


Recombination Component 




1.7 ± 1.3 




eV 


F RRC 

1 cvn 


8.8 ±3 




10- 5 phs-'em- 2 



facilitate comparison 
t lAnders & Grevessd fl989h 

NOTE: Abundance uncertainties are 90% confidence range 



in the spectrum due to the target extent. This broad- 
ening is only slightly wider than the LETGS point- 
source line spread function and can be reasonably ap- 
proximated by a Gaussian profile. In order to test 
the Gaussian spatial broadening approximation, we 
fit for the energy dependence of the smoothing width 
AE cc E a . The best-fit value obtained for a is 2.07 + 
0.15, consistent with spectral broadening of gratings 
in which AA is roughly independent of A, and thus 
AE cc AA/A 2 oc E 2 . 

The best-fit spectral model is plotted over the 
source spectrum in Fig.Q] The most prominent lines in 
the spectrum are those of C, O and Ne K-shell ions. Fe 
L-shell lines, often the best temperature indicators, are 
in fact not clearly identified in the spectrum due to the 
low Fe abundance. The best remaining temperature 
indicators are the line ratios of O VIII to O VII and of 
Ne X to Ne IX, which require here two non-degenerate 
temperature components of kT = 160 eV and kT = 
260 eV. The two temperatures likely represent a con- 
tinuous temperature distribution in the hot bubble. For 
the case of a spherical wind, these range from the 
hottest inner region right behind the reverse shock, to 
the cooler, previously heated, more extended regions 
near the contact discontinuity (CD) that are still much 
hotter than the few-eV outer nebula (see, e.g., Figs. 1 



& 2 in Akashi et al. I l2006h . Much lower temperatures 
(kT < 100 eV) are usually difficult to constrain from 
X-ray spectra, as the emissivities of most bright X-ray 
lines decrease strongly at these low temperatures. This 
problem is augmented here by the significant photo- 
electric absorption (Nh — 2.4 X 10 21 cm 2 ) toward 
BD+30°3639 and by the decreasing LETGS through- 
put with decreasing energy. The present spectrum, 
however, does provide unique temperature diagnostics 
as described in 92.3l and in 33.21 

Relative abundances can be accurately measured 
with X-ray emission line spectra. The measure- 
ment of absolute abundances Az/Ah relative to hy- 
drogen, on the other hand, is prone to systematic un- 
certainties as it requires assumptions for the H con- 
tribution to the bremsstrahlung continuum, which is 
harder to measure than narrow emission lines. In 
the case of BD+30°3639, the low temperature (low 
bremsstrahlung intensity), the low H abundance, and 
the generally non-solar composition of the X-ray 
plasma resulting from core He-burning, all make the 
assessment of Ah even more uncertain. Consequently, 
we will limit our discussion in this paper to relative 
abundances, which are presented in Table Q] The 
present model par ameters can be seen to very well 
agree with those of lYu et al.1 d2009al) . as expected. For 
more details on the uncertainties associated with the 
absolute and re lative abundance determination, see 
Yu et alJ (I2009al) . For the purpose of the present anal- 
ysis, it is sufficient to note the unambiguously high C 
abundance in the hot bubble of BD+30°3639. 

2.2. The Carbon Footprint of BD+30°3639 

The 2-T model described above fits the overall 
spectrum quite well, but fails to reproduce the signif- 
icant emission of what appears to be a dense complex 
of spectral lines in the 25-26 A range. These lines co- 
incide in wavelength with the high-order C VI Lyman 
series. The top panel in Fig. [2] zooms in on the rele- 
vant spectral region showing the data and the best-fit 
2-T model, while highlighting the C VI Lyman series. 
The APEC model includes the C VI Lyman series, but 
only radiative decays from levels with n < 5, n be- 
ing the principal quantum number of the upper level of 
the radiative transition. Thus, we first use the HUL- 



LAC code (IBar-Shalom et al.l 1200 lh to complete the 
spectrum of the Lyman series up to n — 8 explicitly, 
and the decreasing contributions of the remainder of 
the series by extrapolation. The fitted model with the 
complete Lyman series is plotted in the middle panel 



3 




28 29 30 31 32 33 34 
Wavelength [A] 



Fig. 2. — Zoom into the carbon spectral region with 
0.1 A bin resolution. Top: Best-fit 2-T APEC model. 
Middle: 2-T APEC model supplemented with the 
complete C VI Lyman series. Note the failure of 
this model to fit the feature around the C VII RRC 
at 25.30 A. Bottom: The full model (red) including 
emission following recombination with cool electrons 
(thinner brown line). Note that interstellar absorption 
increases drastically with the wavelength across this 
range, which explains the skewed Lyman series line 
ratios. See Fig. 3 for the fluxed model. 

of Fig. |2] Evidently, the agreement of the model with 
the observed data is somewhat improved. 

Even after incorporating the complete Lyman se- 
ries, the most prominent residual feature remains at 
25.30 A. We identify this feature as the narrow radia- 
tive recombination continuum (RRC) of bare C VII 
recombining to form H-like C VI. In the following, 
we consider alternative identifications, but rule all of 
them out. First, this feature need not be confused 
with the N VII Lya line, which is very weak, but 
unam biguously resolve d at 24 .78 A (see Fig. |2}. In- 
deed, IMurashima et al.l d2006l) . using a lower resolu- 
tion spectrum of BD+30°3639, reported a high N/O 
abundance ratio. With the gratings, the two features 
are cle arly distinct and the actual N/O ratio is much 
lower (lYu et al.ll2009al) . Note that an overestimation 
of the N abundance should be expected from any low 
resolution observations that employ a model missing 
high-order C VI lines. Moreover, if a narrow RRC is 
present in the spectrum, it too would be confused at 
low resolution with N VII Ly a and would lead to a 



further overestimate of the N abundance. 

Other potential candidate lines around 25.30 A in- 
clude the emission lines of L-shell ions of mid-Z el- 
ements such as Si, S, Ar, and Ca. The L-shell ions 
of S have no bright emission lines between 25 - 26 A 
( Lepson et al.l l2005) and neither does Si (Lepson et al. 
in preparation). The strongest line of Li-like Ar XVI 
is at 25.02 A dLepson et al. 2003 ). which given the 
LETGS resolving power, is easily distinguished from 
the observed emission feature at 25.30 A. Moreover, 
Ar XVI has equally bright lines at 24.87 and 23.53 A, 
neither of which is observed in the present spectrum. 
No other Ar ion has bright lines between 25 - 26 A. 
Our atomic computations show that Ne-like Ca XI has 
a 2p-4d line at 25.38 A. However, a Ca XI 2p-3d line 
at least five times as strong is predicted at 30.47 A, 
and should be accompanied by the bright 2p-3s lines 
at 35.65, 35.67, and 35.79 A, none of which are ob- 
served in the spectrum. Given the extremely non-solar 
abundances observed in some PNs, we also checked 
the spectrum of Sc. O-like Sc XIV turns out to have a 
relatively bright 2p-3s line at 25.39 A. However again, 
stronger lines of Sc XIV that are not observed are ex- 
pected based on our computations (e.g., a 2p-3d line at 
22.95 A), as are lines of other L-shell Sc charge states. 

Charge exchange (CE) is another process that typi- 
cally enhances high order lines when highly charged 
ions mix with neutral atoms. In the context of as- 
trophysical shocks, unshocked neutral atoms can eas- 
ily penetrate the CD (electro)magnetic barrier and mix 
with the hot gas. If the enhancement is in sufficiently 
highly excited levels of C VI (n > 10), the result- 
ing lines could hypothetically fall below 25.5 A and 
would not be resolved from the C VII RRC in the 
LETGS spectrum. However, there are two strong argu- 
ments against such CE taking place in BD+30°3639. 
First, the UV radiation from the central star main- 
tains high ionization (e.g., O III) in the inner parts 
of the cool nebular gas just outside the CD. Both 
BD+30°3639 and its sister PN NGC 40 feature strat- 
ification in ionization, in which the innermost regions 
of the optical nebula, just outside the hot bubble , 
are the most highly i onized dBrvce & Mellema 1999; 



Sabbadin et al 



2000h . Observations (and an estimate 
of the ionization parameter U ~ 0.03) imply that the 
fraction of neutral species available for CE with the 
X-ray ions is less than 1/1000. The second argument 
against CE in BD+30°3639 is that in order to produce 
an enhanced spectral signature below 25.5 A, the ion- 
ization potential of the neutral atom would need to be 
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Fig. 3. — Fluxed best-fit spectral model for 
BD+30°3639. The small contributions of C VII re- 
combination with cool electrons (as well as the up- 
per limit for the C VI RRC - dashed line) compared 
with the strong line emission originating predomi- 
nantly from hot gas demonstrate how hard it is to de- 
tect recombination from ions crossing CD surfaces in 
astrophysical shocks and how the extremely high C 
abundance of BD+30°3639 (A c /A ~ 40) is crucial 
for the present detection. 



ridiculously low (< 2 eV, see iJanev & Winterl ll985). 
This immediately rules out all of the abundant ele- 
ments in BD+30°3639 including He, C, O, and Ne (as 
well as H) by a large margin and essentially all other 
elements as well. 

An intriguing possibility is CE with dust grains. 
Dust grains can survive UV charging and may have 
low ionization energies (i.e., work fun ctions), which 
could b e of the order of a few e V. Indeed, Matsumoto et al 
d2008l) found evidence of silicate grains in BD+30°3639 
on scales of ~ 4" that coincide with the diameter of 
the X-ray determined CD. However, at the correspond- 
ing distance of 4 x 10 16 cm from the central star, the 
grains inevitably are positively charged by more than a 
few ele mentary charges (e) a s can b e inferred from the 



work of iFeuerbacher et alj (1973). This completely 



reduces the cross section for CE with positive (C VII) 
ions. In short, we deem the CE process as highly un- 
likely to produce the emission around 25.30 A. All of 
this leads us to the conclusion that the 25.30 A emis- 
sion feature must be the RRC of C VII recombining to 
form C VI. 



2.3. The Carbon RRC 

The narrow width of the RRC, which is a direct 
measure of the electron energy distribution in the 
plasma, indicates that recombination occurs with elec- 
trons of a few eV, two orders of magnitude colder than 
those required to ionize carbon to its C VII state. In 
order to measure the exact recombining-electron tem- 
perature, we used the XSPEC redge model, which 
employs the flux density profile 



f E =F R 



(photons s 1 cm 2 eV ') 



kT: 



cool 



(l) 

where F RRC is the total photon flux in the RRC, kT c °° l 
corresponds to the temperature in eV of recombining 
electrons, hv (= Ej+ electron energy) is the photon 
energy, and Ej is the ionization potential of the recom- 
bined ion. The ionization energy of C VI was set at 
Ei = 490.02 eV and the two remaining independent 
parameters were fitted, yielding kT"' ' = 1 .7 ± 1 .3 eV 
and a total RRC flux of F^ c tt = (8.8 ± 3) x 10~ 5 pho- 
tons s^'cm -2 . The upper limit on kT c e °° l is robust as the 
maximum width of the RRC is tightly constrained by 
the data. Conversely, the lower limit is not as well de- 
termined due to the RRC being just broader than the in- 
strumental line spread function expected from the ex- 
tended angular size of BD +30°3639. At th e distance of 
BD+30°3639, d=\.2 kpc dLi et al.ll2002h . we find the 



isotropic rate of C VII recombination events resulting 
in photons of « 25.30 A (i.e., with cool electrons) to 
be 



L cvn 



And 2 F RRC 
r CVII 



(1.5 + 0.5) x 10 



,4(1 



(2) 



For physical consistency, we also include in the 
model the C VI line intensities due to recombination 
onto excited levels (up to n = 8) and ensuing radia- 
tive cascades, with relative contributions that are ap- 
propriate for kT™" 1 . In the recombination model for 
kjcooi -17 e V 5 the unabsorbed relative intensities of 
the RRC, Ly a, Ly p, Ly y, and Ly 5 are 1 .0, 1 .35, 0.22, 
0.08, and 0.04, respectively. These contributions wer e 
calculated using HULLAC dBar-Shalom et alj l200ll) , 
and here added to the model as narrow Gaussians at 
their fixed wavelengths. The final model around the 
relevant spectral region can be seen in the lower panel 
of Fig. [2] including the RRC and the low contribution 
of recombination to the lines compared to the high- 
T emission. Note that photoelectric absorption in this 
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part of the spectrum is significant (A^ = 2.4 x 10 21 
cm" 2 ), and decreases strongly away from the C I edge 
at ~ 43.5 A toward shorter wavelengths, so that low 
order lines of the series (Ly a, Ly [}) are much more 
absorbed than the high order lines. The best-fit flux 
model corrected for interstellar absorption is presented 
in Fig.[3]and the model parameters are listed in TableQ] 

The miniscule flux due to recombination of highly 
ionized atoms with cool electrons compared with the 
strong line emission originating from hot gas makes 
the detection of the former tremendously hard. It is 
obvious from Figs. [2] and [3] that if not for the unusu- 
ally high C abundance (A c /A * 40) in BD+30°3639, 
this detection would not have been possible. Never- 
theless, we did look for more evidence of recombina- 
tion in the spectrum. The RRC of C VI forming C V 
falls at 31.63 A. Due to the low count rate at these 
long wavelengths, only an upper limit can be obtained, 
which is: 3.7 x 10" 4 photons s" 1 cm" 2 (90% confi- 
dence). This flux is high (not very constraining), par- 
ticularly when compared with the C VII RRC flux of 
(8.8 + 3) x 10" 5 photons s 'cnr 2 . The difficulty is due 
to the low effective area of LETGS at this wavelength 
(~ 2 cm 2 ) and the strong absorption near the C I edge. 
Nonetheless, this upper limit still proves useful, as will 
be demonstrated in §|3] The flux in the O VII forbid- 
den line at 22.01 A and the He/3 line at 18.63 A can be 
seen in Fig. Q] to exceed the flux predicted by the hot 
plasma model, which constitutes further, though tenta- 
tive, evidence for recombination as the upper levels of 
these lines are preferentially populated by recombina- 
tion and radiative cascades. The O VIII RRC forming 
O VII falls at 16.78 A, but only an upper limit (90% 
confidence) to the flux of 7.6 x 10" 6 photons s" 1 cm" 2 
can be obtained. Blending with the Fe XVII line at 
that wavelength is possible, but not expected to be sig- 
nificant, since Fe is considerably underabundant. Al- 
though unambiguous detection is possible here only 
for the highly overabundant C, the upper limits are all 
consistent with the coolest (kT a: 100 eV) species in 
the hot bubble recombining with kT » 2 eV electrons 
in the nebular gas. 

The population of intermediate-temperature elec- 
trons (kT > 3 eV) cannot be significant, as it would 
have smeared out the high contrast of the C VII RRC. 
To demonstrate this, we attempted to fit the same 
model as before, but with two C VII RRC components 
instead of one. We note that the quality of the data 
does not warrant a second RRC (see Fig. [2]). Never- 
theless, the main component was fixed at its previous 
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Fig. 4. — Constraints on C VII recombination rates 
from two electron-temperature components as a func- 
tion of kT2 the temperature of the second component. 
The temperature of the first component is held constant 
atkT = 1.7 eV. 



best-fit value of kT\ = 1 .7 e V, while another cool com- 
ponent kT2 was tested for (not to be confused with the 
high temperatures of the 2-T plasma model of 32. 11 1. 
The flux of both components was allowed to vary in 
the fit, as were the other model parameters. In Fig. |U 
we plot lcr flux limits of each of the two RRCs as 
a function of kTi. The 1.7 eV component is signifi- 
cant at all kT2 values, corroborating the existence of 
cool electrons. At high temperatures (kT2 > 20 eV), 
the solution gradually converges to the single RRC re- 
sult. A 10 - 20 eV RRC is possible, although it is 
consistently weaker than the 1.7 eV RRC (Fig. [4]). Be- 
low ~10 eV, the plot shows that there is a high degree 
of confusion between the two components. If we al- 
low both RRC temperatures to vary, a best fit is ob- 
tained for kTi = 1 eV, kT 2 = 10 eV, h = l^j, 
h - 1^0 7 1 40 recombination-photons s" 1 . These do 
not represent distinct components as the temperatures 
are highly confused and interchangeable. Only an up- 
per limit of kT < 20 eV can be obtained for the hotter 
components. The addition of the second RRC (two ex- 
tra free parameters) lowers the reduce d-^ 2 from 0.782 
to 0.775 (Ax 2 = 1 .97), if the method of iGehrell (1986J) 
is used. We conclude that recombination with elec- 
trons hotter than ~ 3 eV is not required by the data (see 
Fig. [2]). However, insignificant recombination with 
electrons of up to ~ 20 eV cannot be ruled out. A 
good probe of the se intermediate tem peratures could 
be O VI emission dGruendl et al .1l2004 . but we are not 
aware of O VI measurements for BD+30°3639. 
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3. Discussion 

In order to produce the observed narrow RRCs, 
highly charged ions have to collide with cool elec- 
trons at temperatures corresponding to only a few eV. 
Narrow RRCs are typical of photoionized X-ray plas- 
mas, where the ionization level is high, but the electron 
temperature remains low; good examples are X-ray bi- 
naries, or Seyfert 2 ionization cones. In contrast to 
such sources, PNs have no hard enough X-ray source 
to ionize carbon up to the C VII state. One might think 
that plasma instability around shocks could entail rapid 
radiative cooling in which ionized atoms recombine 
with cool electrons. However, recombination is one of 
the primary cooling mechanisms, and highly charged 
species such as C VII would have long recombined by 
the ti me the plasma cooled from above 100 eV to a few 
eV dSufherland & Doprtall99"l Fig. 16 therein), which 
rules out instability and rapid cooling as the origin of 
the RRCs. In principle, adiabatic expansion could pro- 
duce RRCs if it is sufficiently rapid that ions do not 
have time to recombine before the plasma cools to a 
few eV, that is, "ionization freezing". In reality, this 
is unlikely. The expansion time depends on the ge- 
ometry, but the confined geometry of the hot bubble 
by the dense nebula does not favor such rapid expan- 
sion. During the expansion, due to the adiabatic re- 
lation T oc « 2 ^ 3 , and since recombination rates (Eq. 



[6]l scale as nT~ 



1/2 



T, recombination would pref- 



erentially occur during the hot, dense phases and not 
during the low-density phases when the plasma cools 
down to a few eV. Hence, it would not produce the 
observed narrow RRCs. Moreover, for the hot gas to 
cool from 100 eV to 1 eV, its density would have to 
decrease by a factor of 1000. Consequently, in or- 
der to produce the observed emission measure of the 
RRC (n^,V ~ 10 53 cm~ 3 ), the size of the source would 
have to be more than two orders of magnitude larger 
than the hot bubble, which can be ruled out. A tran- 
sient ionizing scenario, as expected for shocked gas, 
cannot produce RRCs either, although the initial ion 
temperature would be much higher than the electron 
temperature. This is because electron impact ioniza- 
tion follows electron heating. Even if one invokes a 
yet-unspecified mechanism to ionize the hot atoms (up 
to C VII) before the electrons heat up, the collisions of 
cold electrons with hot ions much more efficiently heat 
the electrons than result in recombination. We con- 
clude that the most plausible physical origin of the ob- 
served RRCs is highly charged species that are heated 
and ionized in the hot bubble (kT > 100 eV) and then 



interact with much cooler electrons (kT ~ a few eV) in 
a distinct medium. 

3.1. Ions Crossing The Contact Discontinuity 

The most obvious reservoir of cool electrons lies 
in the cold nebula, where typical temperatures indeed 
correspond to a few eV. The nebula is photoionized 
by the central UV source, so unbound electrons are 
abundant. The observed recombination can take place, 
therefore, if ions from the hot bubble cross the CD 
plane and penetrate the slow wind where the cool elec- 
trons reside. In the simplified hot bubble picture, the 
outer part of the hot bubble is also its densest and rel- 
atively coolest region (for the spherica lly symmetric 
case see, e.g., 



Akashi et alJl2006i l2007h . This picture 



is consistent with the lowest charge states in the hot 
bubble crossing the CD and recombining, as the spec- 
trum seems to suggest. The narrow RRC and the con- 
sequential limits on intermediate temperatures should 
help constrain the role of heat conduction by electrons 
as well as the role of nebular evaporation into the hot 
bubble. 

Heat conduction by electrons can be supp ressed 
even by magnetic fields as weak as 0.1//G (iSokei 
1994). The viability of such magnetic fields has not 
been established, but if they do exist, they could pre- 
serve the steep temperature gradient across the CD re- 
quired to produce a narrow RRC. The length scales for 
such a gradient are of the order of the electron Lar- 
mor radius, while the Larmor radius for ions can be 
much larger. For particles with mass m and charge 
q at a temperature T, Ri ~ y/mT/q. Thus, ions can 
effectively penetrate the cold nebula across a magnet- 
ically held CD, while the opposite effect of cold elec- 
trons penetrating the hot bubble would be negligible. 
The Larmor radius of a bare C VII ion with kT lw ' « 
100 eV gyrating along a typical interstellar magnetic 
field component of B — 1 fiG (parallel to the CD sur- 
face) is 



/ b r 1 / kT ho ' \ I/2 



R L = 6x 10* 



The corresponding half-circle (cross-back) rotation 
time would be: 



ncM c 

Tr = = 650 

6eB 



(4) 
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where Mq is the atomic mass of carbon. For compar- 
ison, the slow-down time by Coulomb collisions of a 
C VII ion with the ambient cold electrons can be es- 
timated from Equation (5-29) of lSpitzerl dl956t) to be 
(independent of ion temperature) roughly: 



1200 



kT 



cool \ 3/2 



leV I \10 4 cm- 



(5) 



where n c f°i is the number density of cools-gas elec- 
trons responsible for slowing down the fast C ions. 
For the temperature and density of Eq. (0, if the C 
ion is stopped by protons t s ~ 800 s (Spitzer 1956, 
Equation 5-28). Obviously, the strength of the mag- 
netic field and to a lesser extent the density in Eqs. 
Q are fairly uncertain (even to an order of magni- 
tude), but the fact that t s is comparable to tl implies 
that ions crossing the CD have a good chance of be- 
ing stopped by the cool plasma, rather than returning 
to the hot bubble. Chances can be much higher, if the 
magnetic fields are weaker than 1^/G, or absent alto- 
gether. Note that in order for the hot C ions to reach 
the 1.7 eV gas, they need to cross any intermediate re- 
gion (say kT — 10 eV) without being stopped. This 
considerably limits the size of such a region through 
t s (T) and the ion velocity, which for intermediate tem- 
perature kT — 10 eV yields ~ 10 9 cm. 

Compared to these short time scales, the recombi- 
nation of C VII with the cool nebular electrons can take 
several months: 



1 



n e a RR (T e ) 



1.8xl0 7 



kT 



coolX 1 / 2 



leV 



(l0 4 cm- 3 ) 



(6) 



is 



where a RR = 5.5 x 10- 12 cnrV 1 (jfc7'f o 7l e V) 
the radiative recombination rate coefficient. These 
timescales suggest the following plausible scenario: 
ions from the hot bubble (kT hot -100 eV) cross the 
CD into the cool nebula. Many of them slow down be- 
fore they can gyrate back, and thermalize by collisions 
with the nebular electrons, cooling down to a few eV. 
After a much longer time, and perhaps after they dif- 
fuse further upstream, the ions eventually recombine 
with the cool nebular electrons. 

3.2. Temperature Diagnostics at the Contact Dis- 
contiuity 

As seen above, the width of the RRC can provide a 
tight constraint on the electron temperature on the neb- 



ular side of the CD. It would also be useful to measure 
the temperature of the hot plasma close to the CD, but 
just on its hot side. This temperature can be readily 
obtained from the limit on the recombination rate of 
C VI. In a steady state, the rate of recombination with 
cool electrons is balanced by the ion crossing rate. In- 
deed, steady state is reached within the typical recom- 
bination times T rec < year (Eq. |6), which are much 
shorter than the PN age. The measured recombina- 
tion rates I rec , therefore, reflect the corresponding ion 
densities n wn and thus the fractional abundances in the 
hot gas f wn {T hot ) oc «,„„. The C VII to C VI ionic 
density ratio is particularly sensitive to the temperature 
around kT w 100 eV. These densities and temperatures 
refer to the edge of the hot bubble, right at the CD sur- 
face. Taking into account that C VII ions crossing the 
CD will produce two photons, one in the C VII RRC 
and subsequently one in the C VI RRC, one can write 
* fcvn(T hc ") and 7J" oc f CVI {T ho ') + f C vii{T h '"\ 



l cvn 
or: 



fcvn(T ho ') 
fcv,(T ho ') 



rrec 
1 CVII 



rrec _ rrec 
1 CVI 1 CVII 



(7) 



The measured upper limit for 7™. Q2.2t places a 
lower limit on the ratio fcvn(T l "")/fcvi(T lK ") > 
0.31+0.1. The theor etical ionic fraction ra tio as a func- 
tion of temperature (Ma zzotta et al .1998) is plotted in 
Fig. along with the measured lower limit imposed 
by Eq. (0. These results indicate that the temperature 
on the immediate hot side of the CD must be at least 



kT 



hot 



± 4 eV. Together with the accurate mea- 



surement of kT C0 ° l - 1.7 + 1.3 eV from the width 
of the RRC, this implies a temperature jump corre- 
sponding to at least 80 eV across the discontinuity, 
with no significant plasma at intermediate tempera- 
tures. In principle, the O RRCs could have provided 
independent diagnostics, perhaps even a valuable up- 
per limit on T 1 ' ', since the Ac/Ao abundance ratio in 
the hot bubble is well constrained. Unfortunately, the 
present upper limit on the O VIII recombination is not 
restrictive enough to provide meaningful constraints. 

The requirement for a steep temperature jump cor- 
responding to ~ 80 eV across the CD is robust and 
it raises an intriguing question regarding the role of 
magnetic fields and heat conduction. Heat conduction 
would drastically change the temperature and density 
profiles across the CD. Instead of a sharp-contrast in- 
terface, a more gradual conduction front would mani- 
fest the transition between the hot bubble and the neb- 
ula as nebular gas evaporates into the hot bubble. On 
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the face of it, this picture seems to be in contrast with 
the LETGS spectrum of BD+30°3639 and the present 
findings. In addition to the sharp temperature jump 
across the CD implied by the present analysis, the fast- 
wind composition of the hot bubble in BD+30°3639 
argues against a significant role for nebular evapora- 
tion in to the hot bubble, at least in this source dYu et al.1 
2009ah . Note, however, that the penetration of hot 



bubble material into the nebula will alter the chem- 
ical abundances of the immediate inner parts of the 
nebula, eventually changing the composition to that of 
the fast stellar-wind. Such abundance gradients in the 
inner nebular regions, gradually changing from fast- 
wind composition to nebular composition as the dis- 
tan ce from the center in creases, were indeed found 
bv lSabbadinetail (2000) in spatially resolved obser- 
vations ofNGC 40. 



Steffen et al. 



(2008) published numerical simula- 
tions for PN hot bubbles that include heat conduc- 
tion in the magnetic-free limit. Strictly speaking, 
these models cannot be applied to BD+30°3639, since 
they assume solar abundances dominated by hydro- 
gen, while the X-ray gas in BD+30°3639 consider- 
ably departs from solar composition, consists predom- 
inantly of He and C, and is extremely H-deficient. 
Nevertheless , we n ote that the temperature profiles of 
Steff en et al.1 ([2008 ) at the outer edge of the conduction 
front do appear to sustain a sharp temperature drop by 
approximately two orders of magnitude (see their Fig. 
4), perhaps due to the rapid cooling of plasma between 
a few eV and 100 eV. However, the spatial resolution 
of these simulations is of the order of 10 15 cm, while 
the thickness of the transition layer estimated in the 
present work is several times 10 13 cm ( jj3.31 >. Care- 
ful modeling of heat conduction on the small scales 
probed here and under conditions closer to those of 
BD+30°3639 (He and C dominated plasma) will be 
needed before making a more conclusive and quantita- 
tive statement regarding the role of magnetic fields and 
heat conduction in PNs with WC central stars. 

3.3. Further Estimates: Carbon Mass Outflow, 
Densities, and Penetration Depth 

The measured RRC intensity (Eq. |2]l implies a CD 
mass crossing rate of Mcvn = (4.7 + 1.6) X 10~ 9 
M yr _1 , and further dividing by the fractional abun- 
dance of C VII at 100 eV (Yn/„ =0.27. iMazzotta et al.1 
1998b . yields a total recombining C mass rate of Mc ~ 
(1 .8+0.6)x 10~ 8 M yr~' . We recall that approximately 
50% of the hot bubble mass of BD+30°3639 is likely 




Fig. 5. — Ratio of C VII and C VI ionic fractions 
in collisional ionization balance, and as a function of 
the electron temperature (solid line). The horizontal 
dotted lines indicate the lower limit for this ratio of 
0.31 ± 0.1 obtained from the C RRC flux ratio (Eq.[7J). 
The dashed vertical line marks the resulting lower limit 
on the hot bubble temperature in the immediate vicin- 
ity of the CD. 



in th e form of C dYu et alj l2009at iMarcolino et al. 
2007b . If this estimate holds near the CD, it implies a 
mass crossing rate ofM* (3.5 ± 1.2) x 10~ 8 M yr _1 . 
The most recent fast-wind mass outflo w rate estimate 
for BD+30°3639 is 5 x 10~ 7 M yr" 1 dMarcolino et al. 



2007), wh ich is an order of magnitude less than the es- 
timate of lLeuenhagen et all dl996i) . Therefore, the CD 
mass crossing rate is, at the most, less than 10% of the 
mass outflow rate and even a smaller fraction of the 
early PN wind that had an even higher M. To that end, 
the ions crossing the CD seem not to have an apprecia- 
ble effect on the stellar-wind, or the shock dynamics. 

With a simplifying assumption of a smo oth CD ge- 
ometr y that ignores wind clumping (see IPrinja et aD 
2007), the ion density can be inferred from the ob- 
served recombination rate. We assume that ions cross 
from the hot bubble across the CD into the cold neb- 
ula exclusively due to their thermal motion. As in the 
previous section, in a steady state in which the C VII 
crossing rate is balanced by the observed recombina- 



tion rate /; 



era- 



one can write: 



L cvn 



0.5vx5 CD n^ 



(8) 



where v ± = -\/kT hc "/Mc is the average magnitude of 
the C VII velocity component perpendicular to the CD, 
whose surface area is Scd- The factor of 0.5 comes 
from the isotropy of the thermal motion. Eq. © can 
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then be used to express the C VII ionic density on the 
hot side of the CD as follows: 



„hot 
l CVIl 



0.5 



St 



2 x 10 34 cm 



kT 



hot 



100 eV 



-1/2 



(9) 

where we assumed the CD surface area Scd to be 
a perfect sphere with the observed X-ray radius of 
BD+30°3639: R CD = 4 x 10 16 cm (2" at 1.2 kpc). 
Using again the ionic fraction /cv//(100 eV) = 0.27 
yields a C densi ty of n c ~ 2 cm" 3 . Assu ming as before 



and f ollowing Ma rcolino et al.l ([2007) and lYu et al 



(2009a), a non-solar composition of gas depleted of 
hydrogen and dominated by He and C (Ac/Ane ~ 0.4 
by number), we can estimate the electron density in 
the hot bubble to be n e w 3nn e ~ 8«c ~ 16 cm -3 . To 
that extent, the observed recombination rates appear to 
be consistent with the scenario of shocked hot bubble 
gas penetrating the slow wind, and with the theoretical 
models and numerical simulations of such a scenario. 

Of course, the assumption of spherical symmetry in 
these models and in Eq. (0 is probably an oversim- 
plification. Collimated fast winds would imply higher 
densities. Furthermore, the simplified picture that only 
includes crossing of hot ions to the cool side must be 
incomplete as it would breach the plasma neutrality 
and would very quickly induce a strong electric field 
that would inhibit further ion crossings. Hence, there 
must be a reverse-charge current, for example of elec- 
trons following the ions, or nebular ions entering the 
hot bubble. Further discussion of the complex mor- 
phology of the electric and magnetic dynamical struc- 
ture around the CD is clearly beyond the scope of our 
study, but these findings could be useful to constrain 
such investigations in the future. 

Finally, the C VII RRC intensity, along with the 
above estimate of the density, can be used to approxi- 
mate the penetration depth AR of the hot bubble mate- 
rial into the nebula. Using the fact that the recombina- 
tion rate Iryn sca l es with the electron and ion densi- 
ties integrated over the entire emitting volume, which 
in turn can be approximated by a thick sphere (S cd) of 
width AR, just outside the observed X-ray hot bubble 
(the CD), one can write 



A/? 



jrec 
1 CVII 



n e a RR S C D{n C vii) c 

where {ncvn) c °° l is the steady-state average C VII den- 
sity on the cold side of the CD. This average likely rep- 



\cool 



(10) 



resents a density gradient that cannot be constrained 
without a valid model for diffusion and plasma effects 
in the nebula. All we can do at this point is scale AR 



with {ncviiY""' , which can be arbitrarily equated with 
n h °< ofEq.[9] to yield 



AR * 3x10 



13 / n e Y l l kT e 
\10 4 cm- 3 / \leV 



cool \ 1/2 



kT 



hot \!/2 



100 eV, 



cm 



(11) 

At the distance of BD+30°3639, this is only ~ lmas 
on the sky, which clearly cannot be resolved by X- 
ray telescopes even if the actual value of {n cvn ) co ° l 
is much smaller. Such a shell needs to cool rapidly 
enough to radiate away the heat due to the penetra- 
tion of the hot gas. A penetration rate of ~ 10 42 parti- 
cles per second (c.f., Eq. |2), each contributing 100 eV, 
yields a heating rate of ~ 10 32 erg s _1 . If the density 
in the nebula is ~ 10 4 cirT 3 , a shell at a temperature of 

6xl0 47 cm 3 
" x (give n 



kT « 2 eV and with a volume of S cdAR 
would cool at a rate of a t least 6xl0 3 3 
A > 10~ 22 erg cmV 1 , ISutherland & Dopital 11993b 



erg s 



safely above the heating rate. Turning the argument 
around and requiring that the shell be able to expel the 
heat from penetrating ions constrain its thickness to a 
more conservative value of AR > 5 x 10 11 cm. 

There is reason to believe the stellar-wind mate- 
rial continues to diffuse into the nebula much after it 
recombines as spatially resolved images of NGC 40 
show abundance gradients in low-ionization species on 
scales of a few arcseconds, or a few times 10 cm. 
The penetration of the stellar-wind from the hot bub- 
ble deep into the nebula is a different effect from nebu- 
lar evaporation, although both processes similarly defy 
the simplified notion of stark disparity in composition 
between the two media. In contrast with the chemical 
composition, however, the temperature disparity does 
appear to preserve the discontinuous nature of the in- 
terface. 

3.4. Ion Crossing of the Contact Discontinuity 
as a General Phenomenon in Stellar- Wind 
Shocks 

The present discovery of a narrow carbon RRC in 
a PN shock is possible owing not only to the cold, yet 
ionized gas of the ambient nebula, but also to the ex- 
tremely high C abundance in the stellar-wind and the 
hot bubble. It would be natural to suppose, however, 
that hot ions crossing the CD should be a more general 
phenomenon of stellar-wind shocks. Indeed, we found 
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two more cases of narrow RRCs in hot plasma sources: 
both are colliding wind binaries with a C-rich Wolf- 
Rayet (WR) star and an O star companion: y 2 Velorum 
and 8 Muscae. 

y 2 Velorum is a WC8+07.5 stellar binary system, 
which like BD +30°3639 (WC9) i s highly enriched in 
C, O, and Ne. ISchild et all (12004) reported the detec- 
tion of C VII and C VI RRCs in the X-ray spectrum of 
y 2 Velorum, and upper limits to the RRCs of oxygen 
ions. Similar to BD+30°3639, the observed recombin- 
ing gas in y 1 Velorum reflects the WR chemical com- 
position. In the binary case, the O star wind plays the 
role the cold nebula plays in PNs, essentially stopping 
the WR wind and sending back a reverse shock. The 
electron temperatu re of the C VII RRC measured by 
Schild etal]j2004 is 3.3 ±0.7 eV, consistent with that 



of an unshocked O star or WR wind. The RRCs in 
y 1 Velorum persist when the thick WR wind absorbs 
the inner, hot X-ray components implying that they 
originate downstream from the stagnation point along 



the sh o ck front (e.g . , regio n 3 in Fig. 9 of Schil d et al. 



2004). Schild et al. (2004) speculated that the narrow 



RRCs could be due to cold plasma in extended regions 
of the wind, photoionized by hard X-rays from the hot 
wind-collision region. However, they also noted sig- 
nificant problems in finding a model that can account 
for both the high ionization parameter required to pho- 
toionize C to its highest charge state (suggestive of low 
density) and the high emission measure of the recom- 
bining plasma (high density). 



Schild et alj d2004l) also raised the possibility that 



recombination may occur following adiabatic cooling 
of the shocked gas. As discussed in $3] this would 
require the (adiabatic) dynamical times to be consid- 
erably shorter than the recombination times, which is 
highly unli kely in a den s e stell ar-wind. Using the pa- 
rameters of ISchild et al.l d2004l) . we estimate the den- 
sity 1 AU away from the WR star to be n e w 4 x 
10 8 cirT 3 , which implies C VII recombination times of 
T rec * 5000 s (for kT = 100 eV). Adiabatic expan- 
sion times are harder to estimate without knowing the 
geometry, but the dynamical time it takes plasma that 
cools down to kT — 1 eV (C soun d ~ 10 km s -1 ) to 
expand to say ~1 AU is roughly 10 7 s » T rec , As in 
the PN case discussed in §3, if the hot gas expands and 
cools from kT > 100 eV to a few eV, the density would 
have decreased by a factor of ~ 1000. This would re- 
quire an enormous volume to produce the observed 
RRCs. We suggest alternatively that the RRCs in 
y 2 Velorum, in analogy with the case of BD+30°3639, 



are due to highly ionized C atoms that were shocked in 
the WR wind, crossed the CD surface, and interacted 
with unshocked electrons of the much cooler, yet also 
ionized plasma from the O star wind. The high carbon 
abundance of the WR wind, as opposed to the much 
lower abundance in the O star wind, suggests that the 
hot ions must come from the WR side. The wind from 
the O star is also shocked (and hot) close to the stag- 
nation point and therefore this hot-WR-ions/cold-O- 
electro ns interface must exist farther away along the 
bow as lSchild et al.l (12004b inferred from observations. 

identified narrow C 



Suga wara et al 



(120081) identified narrow C VII 
and O VIII RRCs in the X-ray spectrum of another 
colliding-wind binary, 9 Muscae, a WC6+09.5 sys- 
tem fairly similar to y 2 Velorum. They suggested the 
possibility of ions escaping the bow shock layer and 
interacting with the wind, but provided no further dis- 
cussion. Again, we suggest the cool recombining gas 
interacting with the highly ionized plasma in 9 Muscae 
be interpreted as evidence for shock heated plasma in 
the metal-rich WR wind crossing the CD and mixing 
with the cooler O star wind. The bright RRCs in the 
two stellar binary systems and in BD+30°3639 sug- 
gest this kind of mixing and perhaps other microphys- 
ical processes around the CD are common to stellar- 
wind shocks in different systems ranging in size from 
~ 10 13 cm in stellar binaries to 10 16 - 10 17 cm in PNs. 
The high C abundance in these sources obviously fa- 
cilitates the detection of the RRCs of C VII and C VI, 
which are relatively isolated in the spectrum. O RRCs 
are harder to identify unambiguously owing to blend- 
ing with Fe L-shell lines. We suspect, nonetheless, 
that the observed RRCs are traces of shocked plasma 
crossing the CD irrespective of element biases. 

An interesting question is whether similar CD 
crossing takes place in young shell-like SNRs and 
whether it leaves a measurable signature that can give 
an idea of the temperature jump across the CD, as 
it does in BD+30°3639 and in the colliding binary 
winds. The fast stellar-wind colliding with the PN is 
qualitatively analogous to the SN ejecta slamming into 
the circumstellar material in young SNRs. However, 
due to the high velocity of the SN explosion and the 
low density of the circumstellar gas, temperatures on 
both sides of the CD, heated by the forward and reverse 
shocks, are both expected to be high {kT > 100 eV 
and possibly much higher). Additionally, the immedi- 
ate region around the CD represents the oldest shocked 
plasma, so the electrons that were initially cool have 
had time to equilibrate with the hotter protons. Con- 
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sequently, ions crossing the CD would not produce a 
narrow RRC even if they recombine. Moreover, unlike 
in PNs, the influence of strong magnetic fields on the 
CD in SNRs is not in doubt, as indicated by radio syn- 
chrotron emission from thin filaments attributed to the 
CD region. The interstellar magnetic field, and in par- 
ticular its tangential component, can be enhanced by a 
factor of a few at the CD due to field compression b y 
simple charge advection (ICassam-Chenai et al 



or by up to two orde rs of magnitude due to hydrody- 
namical instabilities dJun et al.ll 1 9951) . These fields are 



more than sufficient to sustain the temperature jump 
and perhaps even to totally suppress ion crossing. By 
substituting a magnetic field of 100 fiG and a density 
of 1 cm -3 into Equations. (O and ©, one easily sees 
that a thermal ion gyrating across the CD of a SNR 
will spend far too little time on the cool side of the 
CD (tl ~ seconds) in order for it to significantly cool 
down and remain on the other side (r, ~ days), even 
if cool electrons were present. Indeed, the highest 
quality published X-ray grating spectra (especially be- 
yond 25 A) of young SNRs do not ap pear to feature 



detectable narrow RRC emission (Rasmusse n et al. 



l2Q0lk lBeharetalJl200lh . We therefore conclude that 
narrow RRCs due to ion crossing of the CD are likely 
a unique signature of relatively slow shocks in dense 
stellar-wind sources. 

4. Conclusions 

BD+30°3639 was observed with Chandra LETGS 
for a total of 300 ks providing the first and so far only 
high-resol ution X-ray spe ctrum of a PN. As already 



shown by lYu et al 



( 2009al) . the spectrum can be fit- 
ted by an absorbed, two-temperature plasma model, 
though it likely represents a distribution of tempera- 
tures between ~100 - 300 eV. The spectrum is dom- 
inated by C, Ne, and O emission lines with very lit- 
tle N and Fe. The carbon abundance is particularly 
high, indicating that the X-ray gas primarily originates 
from the present-day WC stellar-wind and making the 
C emission features most suitable for detailed analysis. 

We detect in the spectrum a narrow RRC of bare 
C VII forming C VI by recombination with cool (kT e = 
1.7 + 1.3 eV) electrons. We suspect that a bright RRC 
of H-like C VI forming C V also appears in the spec- 
trum, but the low S/N of that feature allows only an 
upper limit to its flux. We interpret the RRCs as direct 
evidence of penetration of hot, highly ionized plasma 
into the cool nebula past the CD, which separates the 



heated fast wind from the slow wind. This finding re- 
quires a steep temperature gradient, as intermediate 
electron temperatures higher than kT e = 3 eV would 
have broadened the RRC beyond its observed width. 
Using the measured lower limit to the flux ratio of the 
C VII and C VI RRCs, the temperature difference be- 
tween the hot and cold plasma across the interface is 
found to be AkT > 80 eV. Such a steep gradient can 
be sustained by magnetic fields, which would signifi- 
cantly suppress heat conduction. On the other hand, 
magnetic fields may not be needed to preserve the 
sharp temperature drop between the hot bubble and the 
nebula. Nu merical simulations of heat conduction in 
the literature dSteffen et al.l2008 ) do appear to preserve 
a sharp gradient, though their grid is too coarse to re- 
solve the processes discussed above. Simulations with 
the appropriate WC-wind composition and at higher 
resolutions will need to be confronted with the present 
results in order to quantitatively test the role of mag- 
netic fields in PNs. 



0.5) X 
3.5 x 



The measured recombination rate of (1.5 
10 40 s _1 implies a mass crossing rate of - 
10~ 8 M yr , which can be a few percent of the to- 
tal fast-wind mass, implying that the fast-wind gas 
successfully makes its way through the nebula de- 
spite the putative magnetic barrier. The deep pene- 
tration of stellar-wind gas into the optical nebula is 
indirectly supported by abundance gra dients observed 



in the innermost regions of NGC 40 (Sabbadin et al 



2000). This picture is quite different from that of neb- 
ular evaporation into the hot bubble for which we find 
no observable sign. Finally, we note two stellar (WC 
+ O) binary X-ray sources with similar RRC features. 
We take this as evidence of the generality of ions cross- 
ing CD surfaces in slow, dense stellar-wind shocks, re- 
gardless of whether the shock front is on the scales 
of a massive-star binary system or a PN. The hotter 
and more highly magnetized environments of CDs in 
SNRs, on the other hand, are less favorable for ion 
crossing, or for producing narrow RRCs. 
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